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Abstract 

We analize electrical conductivity controlled by hopping of bound spin polarons in disordered 
solids with wide distributions of electron energies and polaron shifts (barriers). By means of 
percolation theory and Monte Carlo simulations we have shown that in such materials at low 
temperatures, when hopping occurs in the vicinity of the Fermi level, a hard polaron gap does 
not manifest itself in the transport properties. This happens because as temperature decreases the 
hopping polaron trades the decreasing electron and polaron barriers for increasing hopping distance. 
As a result, in the absence of the Coulomb correlation effects, in this variable-range variable-barrier 
hopping regime, the electrical resistivity, p, as a function of temperature, T, obeys a non-activation 
law: In (p/po) = (t/T\ with p = 2/(d + 2), where d is the dimensionality of the system. It differs 
from the standard Mott law for which p = l/(d + 1). Also, we studied the effects of upper and 
lower boundaries in the polaron shift distribution on hopping conduction, which may result in a 
partial re-entrance of the hard polaron gap. We discuss possible applications to the problem of 
giant negative magnetoresistance in dilute magnetic semiconductors and nano composites where for 
paramanetic materials p = 3/(e? + 2). 

PACS numbers: 72.20.Ee, 75.30.Vn, 75.50.Pp 



I. INTRODUCTION 



Electrical conductivity due to variable-range hopping (VRH) in doped crystalline and 
amorphous semiconductors has been a subject of intensive experimental and theoretical 
studies for more than three decades 0, [| [|. The temperature dependence of the VRH 
conductivity is described by stretched exponents a oc exp (—T/Tq) p and depending on the 
value of p can be attributed to one of the two major regimes such as the Mott (p = 1/4) 
and the Shklovskii-Efros (p = 1/2) regimes [[I], Several observations of giant and colossal 
negative VRH magnetoresistance in systems with strong spin disorder such as dilute mag- 
netic semiconductors (DMS) have been reported recently R [| [TJ. It seems rather natural to 
attribute VRH magnetoresistance to spin polarons that are believed to be a major cause of 
negative magnetoresistance in DMS || |9|, [LOU . However, a theory of VRH magnetoresistance 
has never been addressed in the literature, partially, because of the belief that spin-polaron 
hopping may lead only to a simple activation dependence of the conductivity which mimics 
nearest-neighbor hopping fi~l| . We will demonstrate that in the systems with significant 
disorder of the magnetic energies spin-polaron hopping is described by a stretched exponent 
similar to that of the Mott behavior. 

The study of electrical conductivity in solids, which is controlled by polaron hopping, 
has centered so far on materials where the polaronic part of the activation energy is fixed 
|1|. [|, [R], [U], [12], [13|, [L4|]. Indeed, the lowering of the free energy of the system due to atomic 



spin alignment |9j ( or lattice polarization |13| ) in the vicinity of the trapped carrier, the 
so called polaron shift, and its associated hopping barrier are generally determined by the 
spatial extent of the localized states between which the hopping occurs. This is usually true 
no matter how strong the coupling between a localized carrier and magnons (or phonons) 
and how many of them are emitted or absorbed during the hopping event |]15| . (Some- 
times, non-linear polaron effects may lead to an additional self-localization of the bound 
state resulting in the actual collapse of its wave function |§ . These situations are out of 
the scope of the present paper for they would hardly be favorable for the hopping process.) 
That is why in standard lightly doped semiconductors the polaron hopping barrier can be 
considered fixed 0. (This is also correct in heavily doped and compensated nonmagnetic 
semiconductors, where at low temperatures conductivity is controlled by electron hopping 
between aggregations of impurities of a fixed size close to the screening length H.) In all 
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these cases, the bound polaron spectrum has a hard gap resulting in a predominantly acti- 
vation temperature dependence of conductivity even in the variable-range hopping (VRH) 
regime 
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Recent technological developments in producing novel dilute magnetic semiconductors, 
especially GaAs : Mn, and magnetic nanostructures, such as GaAs/ ErAs, made it possible 
to observe electrical conductivity due to hopping via magnetic impurities and nano-islands 
|16| thus revealing a rich interplay of magnetic and electrical transport properties. For in- 
stance, Mn atoms play a dual role as acceptor sites and magnetic ions in GaAs resulting in 
the very unusual magnetic properties of this material fg, |17]. Giant negative magnetoresis- 
tance has been observed [||, || [7| below the metal-insulator transition point of GaAs : Mn, 
which clearly displays the VRH behavior with no sign of the hard magneto-polaron gap. The 
absence of this gap in the electrical conduction experiments can be understood by taking 
into account local fluctuations of the bound polaron shifts. These, in turn, may be related 
to the distribution of exchange couplings between Mn ions and holes as was suggested in 



Ref. [17|. In addition, it is quite possible that the fluctuations of the local spin polaron 
shifts are not correlated with those of the "bare" electron energies. In heavily doped and 
compensated samples [0] of GaAs : Mn, the former can be determined by the total number 
of magnetic impurities in the localization volume while the latter may depend on the elec- 
trically active ones. For ErAs nanoislands embedded in GaAs matrix [T(| , the distribution 
of polaron shifts follows the distribution of the islands' volumes. Meanwhile, the "bare" 
electron spectrum of the islands is apparently determined by the distribution of the islands 
sizes and shapes. We will show that fluctuations of the local (lattice and spin) polaron 
shifts in certain circumstances may wash out the hard polaron gap, leading to a substantial 
modification of the standard Mott law for the electrical conductivity p(| even in the absence 
of Coulomb correlation effects. 

Bound spin (magnetic) polarons, in contrast to their lattice counterparts, possess polaron 
shifts that, as well as their distribution, can be tuned by applied magnetic field thus resulting 



in giant or even colossal magnetoresistance [[10], |ll], [14], [16|] . Externally controlled switching 
off and on of the magneto-polaron effect also allows one to better identify the hopping 
mechanism in magnetic materials. 

The main goal of the present paper is to develop a consistent theoretical approach to the 
description of variable-range hopping of polarons in materials where disorder gives rise not 
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only to the electron-energy distribution but, also, to the distribution of local polaron shifts. 
This approach is based on percolation theory and Monte Carlo simulations. In Section [TT], we 
will present simple qualitative arguments concerning the main results of the present paper. 



the Monte Carlo simulations of the polaron hopping problem aimed at finding the universal 
constants and functions that appear in the applications of the percolation theory. Effects 
of the low boundary in the distribution function of the polaron barriers will be considered 



semiconductors and nanostructures, a phenomenon that in some materials can be controlled 
by variable-range variable barrier hopping of the spin polarons. 

II. QUALITATIVE CONSIDERATION 

We will start from a qualitative description of the low-temperature electrical resistivity in 
a system where localized states have electronic energies randomly distributed in the vicinity 
of the Fermi energy Ep. The derivation of the original Mott's law is based on the assumption 
that the density of states in such a system near Ep is constant [l], pOfl . In addition to 
this assumption, we will suppose that the states in questions possess randomly distributed 
polaron barriers as well. Our goal is to find out how this factor would modify the standard 
Mott's law for DC hopping conductivity. 

In the original Mott's derivation, it is assumed that at a given temperature T most of 
the hopping events involve random sites with one-electron energies belonging to a strip of 
width e around the Fermi level. In our case, we have to account for an additional random 
distribution of the polaron barriers. Thus we generalize the Mott conjecture and introduce 
a hyperstrip of width e and height W, where W is a scatter of the polaron barriers. The 
probabilty of a typical hopping event for two sites belonging to this hyperstrip reads: 



where a is the localization radius of the electron wave function, T is measured in the energy 
units, and r is the average distance between the sites of the hyperstrip: 



A rigorous theoretical derivation of these results is given in Section |TJ. Section [TV] deals with 



in Section 0. In Section VI we briefly discuss giant negative magnetoresistance in magnetic 




(1) 



r ~ N~ 1/d = (GeWy 1/d . 



(2) 
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Here N is the concentration of the sites in the hyperstrip, d =2 or 3 is the spatial dimension 
of the system and G is a joint density of (polaronic) states (DOS), such that G (eo, Wo) dedW 
is a probability, per unit volume, to find a site with e in the interval (e , e + de) and W in 
the interval (Wo, W$ + dW). (As in the original Mott's derivation jj], |S{J, we assume that 
G is constant. However, this is not true if the Coulomb correlations are taken into account 

00 

After substituting (g) in flU), it is easy to find that the probability ([!]) as a function of e 
and H 7 will reach its sharp maximum at 

/ 2T \ rf /( d+2 ) \2/(d+2) 

£ - = ^ 2 =(^) =r(7f/T) . (3) 

The hopping conductivity, a, is proportional to w(e opt , W opt ). Therefore, substituting the 
optimal values (Q) and (0) in (fj) yields the following expression for the hopping resistivity: 



p= 1/a = p exp 



T Q [d) /T) 



2/(d+2) 



(4) 



where 



and /3q is a numerical factor of the order of unity. Expression (^) constitutes the main 
result of the present paper. (In the next section, we will obtain this result by means of 
a more rigorous percolation theory approach. Also, by using the Monte Carlo simulation 
described in Section 4 we will be able to find the numerical factor f3^ in Eq. (H).) It can 
be seen that, in accordance with the above derived formula, in the 3D (2D) case the well 
known Mott's exponent 1/4 (1/3) (see 0, [2tJ) should be substituted with a somewhat larger 
exponent of 2/5 (1/2). 

As the temperature goes up both the width and the height, (|3|), of the optimal hyperstrip 
increase while the optimal hopping distance @ decreases. It means that with increasing 
temperature the polaron chooses to jump over increasingly smaller distances at the cost 
of overcoming the increasingly higher and higher polaron (and electron energy) barriers. 
That is why we call this regime a variable-range, variable-barrier hopping (VRVBH). In this 
regime, in contrast to the nearest-neighbor hopping (NNH) of polarons, both the lattice 
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JH, |3| or the magnetic ones QTOJ , the polaron hard gap is washed out by the fluctuations of 
the polaron shifts (barriers). 

If, however, there is a maximal polaron shift in a system, W may[ , then at high enough 
temperatures, such that T > T x , where 



n 



(d) 



t {d) w 



(6) 



the optimal polaron shift @ becomes than VFmax. (Here t$ is a numerical factor of the 
order of unity, which will be determined later by Monte Carlo simulation.) In this case, 
the combined DOS, G, should be substituted with an ordinary electron DOS, g, and the 
VRVBH law (|) with a slightly modified Mott's law for polaron VRH 0, [| 



p = po exp 



, „ x l/(<2+l) 

T$/T) + r d W max /T 



(7) 



where the second (hard-gap activation) term under the exponent sign is small compared to 
the first (Mott's) term. Here 



i(d) _ Pm 



ga a 



- M 



(8) 



where and (3$ are numerical factors of the order of unity, which will be established later. 

At low temperatures, both the electron, g, and the combined, G, DOS cannot anymore 
be treated as constant in the vicinity of the Fermi level due to the appearance of the soft 
Coulomb gap in the electron spectrum 0]. The latter is related to the long-range Coulomb 
correlation effects in the spatial distribution of electrons above and holes below the Fermi 
energy. The maximal width of this gap can be estimated as [H 



i/(d-i) 



(e 2 / Ka) d /T 



(d) 



i/(d-i) 



(9) 



where k is the dielectric constant of the material and e is the electron charge. This gap will 
modify the VRVBH law (|J) at low temperatures when its width exceeds the width (|[) of 
the optimal strip: ec > e op t, i-e. when T < T% with 



if ~ e c (e c /n a >) 



n(d)\ 



2/d 



(10) 
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where Tq is given by expression (|5|). It can be seen that the temperature interval (T 2 < 
T < T[ d) ), where the VRVBH regime (W) exists, may be sufficiently large if the polaronic 
effect is strong enough. Based on the evaluations of G, presented in Section 6, it is possible 
to show that the required condition << is equivalent to W max >> ec- 

III. CALCULATION OF RESISTIVITY. A PERCOLATION THEORY AP- 
PROACH 

Let us consider two sites (quantum dots, impurity centers or their aggregations, etc.) i 
and j separated by a distance r^. To find the bound polaron hopping rate in the so called 
nonadiabatic limit J12], [13] one can employ the semiclassical approach developed in [12|, [T3| 
for the lattice small polarons (SP) or in for the bound magnetic polarons (BMP). In 
both cases, however, is possible to show (see Appendix) that if i — > j tunneling transition is 
governed either by strong multiphonon coupling (for the lattice SP) or by thermodynamic 
fluctuations of the local magnetization (for the BMP), when Wi + Wj ^ |e,- — efj| , the hopping 
rate 



jij = 7°- exp 



2r i3 (Wt + Wj+Ij-Itf 



(11) 



a AiWi + W^T 

where a is the spatial-decay length of the localized wave function into the semiconductor 
matrix, q = e/ (0) — W\ is an equilibrium energy of the /-th site with the localized carrier; 
€i (0) is the "bare" electron energy of the carrier-free state and Wi is the bound polaron 
shift. The latter describes the lowering of the free energy of the system due to either lattice 
polarization or to atomic spin alignment in the vicinity of the trapped carrier. (Usually 
Wi is inversely proportional to the number of polar or magnetic atoms in the localization 
volume 0,0.) 

If, however, electron-phonon (electron-magnon) coupling is not strong, i.e. if Wi + Wj < 
1 7j -£i\ , the single-acoustic-phonon assisted tunneling of Miller- Abrahams type || prevails 



over the resonant one described by Eq. (JTTJ) . Then (see 0, (TO]] for details) 



7y oc (N P + 1/2) ± 1/2, (12) 

where 
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-l 



(13) 



is Planck's distribution function. Here, in Eq. (|TTJ) , "+" and "-"correspond to the cases of 
phonon emission (e,- < and absorption (e} > e*), respectively. 
The effective impedance between the localized sites in question @ 



% = TT\ T\ (14) 



T 

e 2 7y/i (1 - fj) 

where fi — {1 + exp [(?; — _Ep) /T]} _1 is Fermi's distribution function of the bound polarons 
with Fermi level Ep. From ( |TTD - fli~4] ) it follows that 



• y-y, xexp(^ • ^-7-^)- (15) 



where 



and 



o, r,-r, >w i + Wj 

(Wi + - ^ - q|) 2 /4 (Wi + Wj) ,\e i -e j \<W i + W 3 
It should be mentioned that Eqs ([15]) - (0), though resembling those previously obtained 
in |T0| , pTT]] , differ from them in one important point. Here the modified BMP energies 
?; = ti (0) — W\ play role of the "bare" electron energies e\ (0) used in pi 0. In the absence 



of fluctuations of polaron shifts (Wi = Wj = W), this substitution is irrelevant because the 
same constant polaron shift W is incorporated into Ep [[|. It is, however, substantial if 
these fluctuations play a significant role as they do in the case of variable-range hopping of 
the bound polarons. 

Now we are in a position to calculate the electrical resistivity of the system by using the 
percolation theory approach 0. We start from the following bonding criterion (see Eq. ([15]) 

C ; ^ • < C (18) 

Let us introduce the maximum values of e', W, and r, compatible with this criterion, 



S 



e opt = T^ W opt = 4T£, rv = </2, (19) 
as well as the following dimensionless variables: 



r g -jgg W l / Qn N 

sj=— , * = > W ' = U^ (20) 

' opi c op i •'►'opt 

In terms of these variables the bonding criterion flOf) can be rewritten as 

+ + (<p) < 1, (21) 

where 



<% = 2-(N + N + ^-^l), ( 23 ) 

and 



M „ ) = < °. P 1* - «,i > 4 («,+«,)_ (24) 

[4 (wi + wj) -ip\6i- 5j\} 2 /16y? (wi + Wj) , <p\5i- Sj\ <A{w i J r Wj) 

Here a dimensionless parameter <p accounts for a possible existance of an upper boundary, 
Wmaxj in the distribution of the polaron shifts W\. If, however, the optimal polaron shift 
W opt = 4T£ < Wmax, this boundary does not play any role, and one should put ip = 1 in 
(f2~4j). This is true at low temperatures, T < T± @, otherwise (p > 1. 

Now the following percolation problem can be formulated in the d+2 dimensional space. 
There are randomly distributed sites in this space with the dimensionless concentration 

4 d) (0 = G ■ 2l opt ■ W opt ■ r d opt = 2 3 ~ d Ga d TY +2 i (25) 

such that each of these sites, I, has two random parameters, 5i and wi, uniformly distributed 

inside the hyperstrip with \Si\ < 1 and < wi < 1. Given the bonding criterion (gTj) 

with ip = 1, one should numerically find a minimal (threshold) dimensionless concentration, 
„(<*) _ „(<*) ftM 



[Cio'ji a ^ which the percolation onset takes place. If such a concentration is 



found then, in accordance with Eq. fl25|), the critical exponent 



Ad) 

ScO 



2 3 - d Ga d T 2 



l/(d+2) 



(26) 



and therefore the resistivity 



p = Po exp(eS ) ) (27) 
is given by the previously obtained Eqs (H) and (Bj), where the dimensionless factor 

} = 2 d - 3 ^ } - (28) 

If there is a maximal polaron shift W max then at high temperatures T > (0), such 
that = 4T^q' ) > VFmax, all the polaron hopping barriers are accessible. In this case, 
GW opt = g and the dimensionless concentration (^5|) should be substituted by 



(0=9- ^o P t ■ r d opt = 2 1 - d ga d Te +1 . (29) 

Then given the bonding criterion ( |21| ) one should find the threshold critical concentration 
^cf (v 9 ) = vf^ (£i?J m the (<i+l)-dimension space for any dimensionless temperature if = 



4T^/W m3X > 1. 

As far as this concentration is found, the critical exponent can be expressed as 



2^ (g) 
ga d T 



i/(d+i) 



(30) 



So, the parameter <£> > 1 can be found as a solution of the following transcendental equation: 



(f = t 



l/(d+l) 



(31) 



Here a dimensionless parameter 



4 T d T ; 



l/(d+l) 



> t 



(d) 



(32) 



has been introduced where the Mott's parameter is given by Eq @ with 



= 2^? (oo) . 



(33) 
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(The minimal value of the parameter fl32|) , £q , which defines the transition temperature 
@, can be found from Eq (|3l|) by setting there </? = 1.) 
For the high-temperature VRH regime (T > ), the resistivity can be presented in the 
following form 



?=A,exp(^)=,„e X p[(T<?/T) 1/W+1) /<«>(*) 
It obeys a modified VRH Mott law with 



(34) 



f[ d) (t) = <p (t) ft (35) 

being a universal function of the dimensionless parameter (^), where the function (p(t) 
itself is a solution of the transcendental equation (0). Function tp (t) can be found given 
the threshold concentration, function of </? > 1. In particular, from ( pT|) - ( ^4[) it 



follows that due to the fact that Ay- > and is a dropping function, z/^ (<£>) / i/^ (oo) > 1 
and is a dropping function too. As a result, f± (t) > 1 is also a dropping function. At high 
temperatures, such that T >> T\, parameter t » 1. Then the solution of Eq. fl3"T|) can 
be presented as (p (t) ~ t and therefore /i (£) ~ 1 in this limit. In this case, the resistivity 
is governed by the conventional Mott's law. (Also, from Eq fl24|) it directly follows that 
Ay — > when (p — > oo. In this limit, the bonding criterion (|21| ) coincides with that of the 
standard VRH problem [[|.) 

It should be noted that the condition T = T^(0) is equivalent to £^jp = i.e. the 
transition from Eq. (^) to Eq. fl3"4"|) takes place when the widths of the corresponding 
optimal energy strips, ? op t are equal to each other and, simultaneously, the optimal polaron 
shift W opt = W max . 



IV. POLARON HOPPING. MONTE CARLO SIMULATIONS 

A Monte Carlo simulation procedure has been developed and implemented, which aimed 
at finding numerical factors /3q , fyf^ and r d in Eqs ([|) and (|7|), as well as the universal 
function f± (t) and parameter in Eqs ([£|) and (^) . This procedure can be described 
as follows. In a cube or square of a size L, N randomly distributed sites are planted with 
an average concentration vl = N/L d , less than the critical one. Each of these sites, I, is 
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ascribed two random parameters, Si and wi, such that \Si\ < 1 and < wi < 1. Any two sites, 
i and j, are considered connected (bonded) if they satisfy the pairwise bonding criterion ( |24|) 
with if > 1, where is the distance between these sites measured in the same units as L. 
Two sites are said to belong to the same cluster if there is a sequence of connections from 
the first site, through other sites, to the second one. (Each cluster has its own identification 
number.) To test for percolation, boundary layers of thickness b are introduced. |2], IE \. For 



a given if > 1, the code checks if a percolation cluster has been established, which has at 
least one site belonging to each of a pair of opposite boundary layers. If it has not, new 
sites are added in small increments and the clusters updated until percolation is detected 
at a critical (threshold) concentration v c l ((f). Obviously, the latter depends on a given 
realization of the percolation cluster, so the procedure was repeated from 50 to 100 times in 
order to find an average critical concentration (i> c l (f))- 

Thus found critical concentrations for several increasing sizes L have been used to de- 
termine a thermodynamic limit (L — > oo) of the critical concentration, v c (if) , by using the 
following power-law extrapolation: 

(v cL (f))=v c (f) + AL- a , (36) 

where A and a are constants. (It has been found that the distribution of v c l (if) is Gaussian 
with dispersion that tends to zero as the size of the system increases.) 

Fig. |I] clearly demonstrates that, for example, for (f — 1 and f> = 1000, the average 
critical concentrations as functions of the size L converge to their corresponding thermody- 
namic limits, v c (if) , that do not depend on the thickness of the boundary layer b. (In our 
simulations, the latter has been chosen to be close to 0.2z/ c (if) and 0.1z/ c (if) for any 
given if > 1. And, as has been explained in 0, for any given L, the thicker boundary layer 
corresponds to a systematically smaller critical concentration.) 

In particular, we have found that the critical concentration vf^ (1) = that determines 
the critical exponent (26) describing the conductivity @) in the VBVRH regime, is equal to 



42.19 ±0.25 and 31.25 ±0.18 for 2D and 3D systems, respectively. This yields the following 
values of the numerical factor (^), which specify parameter T in Eq (|]): (3^ = 21.09±0.12 
and 4 3) = 31.25 ±0.18. 

Also, the above described procedure has allowed us to establish values of the critical 

(2) (3) 

concentrations in the limit of large if. (oo) = 6.88 ± 0.08 and v K cl ' (oo) = 4.24 ± 0.03. 
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FIG. 1: Average critical concentrations, (v c l (<£>)), as functions of size, L, of the 2D and 3D system 
calculated for two different thicknesses, b, of the boundary layer for tp = 1 (the upper four curves) 
and for cp = 1000 (the lower four curves). The fitted curves are obtained by means of Eq. (p£ 



These, in turn, determine the numerical factors (|33|) in the Mott's law parameters ([|): 
P$ = 13.76 ±0.16 and /3$ = 16.96 ±0.12. They can be compared with previously obtained 
values = 13.8 ±1.2 and &g = 21.2 ±1.2 § (Our somewhat smaller, and more accurate, 
value of /3^,can be justified by using up to 3 x 10 5 sites in the present simulations vs the 
previously used 1.5 x 10 2 sites 0.) It should be mentioned that our values of are 
in a surprisingly good agreement with the ones {(3 M = 13.3 and r)J = 16.4) obtained by 
Ioselevich [[HJ who applied an approximate multicomponent percolation criterion to analysis 



of both the NNH and VRH problems. 

Fig. |] shows the calculated universal function (|35]) which describes a continuous transition 
from the VRVBH regime of polaron conductivity (|J) to the standard Mott regime. For t > to, 
it can be approximated as 
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FIG. 2: Universal function (t) found by means of the Monte Carlo simulations (circles). The 
solid curves correspond to extrapolation formula p7| ) 

TABLE I: Parameters of the universal function (t). 





a 
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c 


to 


00 


Pm 


T 


d = 2 


0.967 


0.756 


-0.157 


0.546 


21.1 


13.8 


0.189 


d = 3 


0.982 


0.694 


-0.178 


0.607 


31.2 


17.0 


0.174 



(t) = a + br 1 + ct~ 2 . (37) 

Fitting parameters for 2D and 3D well as parameters = 6/4 from Eq. (^) and 

other relevant parameters can be found in Table |I[ 
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V. LOW TEMPERATURES. MINIMUM POLARON BARRIER EFFECT 



In certain disordered systems, such as doped semiconductors close to the metal-insulator 
transition, the ditribution of polaron shifts, being broad on its high side, may be very sharp 



on its low side [IT|. It means that a system may have some minimal polaron shift, 1^. 
In this case, the law (Q) is also violated at low temperatures when the optimal polaron shift 
becomes smaller than the minimal one: W opt 



4T£$ < W min , i.e. at T < T 3 , where 



T^^(W min /W m ^ d+1)/d Tl d) . (38) 

where is given by Eq (The regime described by Eq. (^) can be observed within a 
temperature interval maa^T^, T^) < T < that clearly exists if (W max — W min ) fee >> 
1.) At such low temperatures, W m \ n /2 is the only polaron barrier accessible and the "bare", 
ei (0) or the modified, = t\ (0) — Wmin, electron energy is the only fluctuating variable. 
The connectivity criterion ( Jl8| ) can be rewritten as £y- = + A£y < £ + A£y, where 
= W min /2T and 



(39) 



Here 



A 



(e;-e~)78W n 



/2, 



e}| > 2W n 



\e l -e j \<2W n 



(40) 



After introducing the maximum values of ei and ~r*i (see (|T9|), compatible with the criterion 
& < £, as well as the dimensionless variables if / and Si (Eq. (p0|) ), the following percola- 



tion problem can be formulated (see also [] 1|). For any given dimensionless temperature, 
6 = 2T^2 /Wioia, one should find the critical dimensionless concentration, u& (9), when 
percolation first appears given the connectivity criterion 



% + + A° (9) < 1 



where 



(41) 



-1/0, 

(5; - <^) 2 /16 - |5i-5j|/2, 



>4 
0|^-<5,-| < 4 



(42) 
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and assuming that Si are uniformly distributed within the strip \Si\ < 1. When such critical 



concentration is known, then the critical exponent, £^?> * s gi ven by Eq. fl30| ) where (y?) 
should be substituted by v& (6) . The resistivity can be presented in the following form 



id) 



p = p exp 



where 



IT 



2 T d T 



■drp{d) 
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i/(d+i) 
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IT 



(43) 



(44) 



Here T^f is given by Eq. (|§|) and (tj) = 1 + 6 (77) . Function 6 (77) is a solution of the 
transcendental equation 



where v& (6 



oo 
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e = ri 



V. 



% (0) 



id) 
J c2 



OO 



l/(d+l) 



(45) 



(if — > oo) has been calculated in the previous section. 



As it follows from Eq. (pffi), A°- (9) < and is a rising function. As a result, 



v. 



(d) 
c2 



(0) / i/£ (oo) < 1 and is a rising function. Therefore, in Eq. 



/i (^) i s a rising 



function as well. In addition, if 6 — > 0, z^2 ($) I u c2 (°°) ~~ * S<f < 1- Then at low tempera- 



(d) 



tures , when 1, / 2 (d) (77) ~ 1 + 



1] which yields |TT 



p = po exp 



2T 



+ 



^ A/ 

T 



i/(d+i)' 



(46) 



This almost activation type of conductivity (for the second term in the exponent is much 
smaller than the first one) will be observed if T << T 3 . 

The Monte Carlo simulations allowed us to calculate the universal function (77) in Eq. 
(531), which can be approximated as 



f!f ] (77) = 1 + 77 [1 - A exp (-B77)] . (47) 

Fitting parameters for this formula and constant <^ = (1 — A) d+1 from Eq. ( |46| ) are presented 
in Table [TJ. 

And finally, if — W min is smaller than the width e opt (^j of the Mott strip, the 

VRVBH regime, described by Eq. (^) ceases to exist. The fluctuations of polaronic shifts 
are not important anymore (W max ~ W^mm), and, if the Coulomb correlation effects are 
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TABLE II: Parameters of the universal function (77). 



A 



B 



d = 2 



0.137 



0.0759 



0.643 



d = 3 



0.109 



0.0612 



0.630 



not substantial (W max >> ec), the resistivity is described by Eq. ( |43| ) at all relevant 
temperatures. It is interesting to mention that in the high-temperature case (t ~ t] » 1,), 
when /j (£) ~ 1 and (77) — rj, both formulas (|34D and fl4~3|) lead to practically the same 
modified Mott law. 

VI. APPLICATION TO GIANT MAGNETORESISTANCE 

In the previous sections, we described a new mechanism of polaron conductivity controlled 
by variable-range hopping which leads to a non-activation temperature dependence 



with p = 2/(d + 2). This result, which differs from the standard Mott's p = l/(d + 1), 
can be justified for both the lattice and the spin polarons if their polaron shifts (barriers) 
fluctuate independently from their bare electron energies. However, the case of bound spin 
(or magnetic) polarons is of special interest because not only the polaron shifts but their 
distribution can be altered under the action of the external magnetic field, H. It can be 
shown, that, indeed, the combined density of polaron states, G, which defines parameter T in 
Eq. (f4"gp (see Eqs (Q) and (|5|)), will decrease with magnetic field. As a result, the resistivity 
will sharply (exponentially) decrease with magnetic field, so that Eq. (f|) also describes 
giant (or colossal) negative magnetoresistance observed in dilute magnetic semiconductors 
and nanostructures. 

To demonstrate this, let us suppose that, for the sake of simplicity, the "bare" electron 
energies e and the polaron shifts W and are not correlated and both are distributed uniformly, 
each in the interval (e — Ae, e + Ae) and (W^in, W mauX ), correspondingly. Then the density 
of the "bare" electron states 




(48) 



17 



.0, |e- el ^ Ae 

9=< , > ( 49 ) 

fl'o = N t /2Ae, |e-e|<Ae 

while the density of the BMP shifts 

, 0, \W-W\^AW 
9w={ ' , ' , (50) 

N t /2AW, \W-W\< AW 

where N t is the concentration of the sites and 2 AW = W min — W max , W = (W m in + W m3X ) /2. 
For the chosen distributions, 

G ~ — — ^° — 9o9w 

~ 2Ae ■ 2AW ~ 2AW N t 1 ' 

given A?~ Ae >> AW. 

For bound spin (magnetic) polarons, their shifts, W, proportional to the magnetic sus- 
ceptibility, x (H, T) [g, 0, fTf. The same is true for the width, 2 AW = W min - W max , of the 
shift distribution fl5Dp. Therefore, by aligning the atomic spins the applied magnetic field 
not only reduces the polaron hopping barriers but also makes their distribution sharper. 
Then, from Eqs. (|5|) and (^TJ) it follows that parameter Tg^in expression ([|) 

ff(tf)oc VxjHTr). (52) 

It substantially decreases with magnetic field when the magnetization approaches its satu- 
ration value thus leading to giant or colossal negative magnetoresistance. 

If the material is in the paramagnetic phase and if the anti-ferromagnetic coupling be- 
tween the magnetic atoms is suppressed, then the susceptibility in Eq. ( |5^ ) obeys the Curie 
law: x {H, T) oc 1/T. In this case, the exponent p = 2/{d + 2) in expression f [48|) describing 
VRVBH regime should be substituted with a somewhat larger quantity p — 3/(d + 2). It 
can be seen that in this regime for 3D systems, 0.4 < p < 0.6. That kind of dependence 
can be easily confused with the standard p = 0.5 typical for electron VRH in the presence 
of Coulomb correlations ||. (The effect of these correlations on VRH of the spin and lat- 
tice polarons will be considered elsewhere. In particular, it can be shown that in this case 
p — (d + 1) / (2d + 1) or p — (d + 2) / (2d + 1) in the paramagnetic phase.) However, at high 
magnetic fields (for spin polarons) the p = 2/(d + 2) or p = 3/(d + 2) VRVBH dependence 
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(for the spin polarons) will be replaced by the p = l/(d + 1) Mott's law while at the same 
conditions the presence of the soft Coulomb gap will result in the p = 0.5 dependence. 

Giant negative magnetoresistance in the variable-range hopping region apparently has 
been observed in recently discovered dilute magnetic semiconductors, such as GaAs:Mn 
|5|, H 0] and MnGe |f2TJ below the metal-insulator transition point. In all these materi- 



als, magnetoresistance is described by Eq. with p = 0.25 4- 0.56 and T significantly 

decreasing with magnetic field. Additional, more detailed measurements of the transport 
and magnetic properties in question are needed to more specifically identify and quanti- 
tatively describe the hopping mechanism in these materials. Also, the effect of the soft 
Coulomb gap should be taken into account when describing electrical conductivity governed 
by variable-range polaron hopping at very low temperatures. 

VII. CONCLUSIONS 

We have examined variable-range hopping of polarons in polar and/or magnetic disor- 
dered materials. In addition to electron impurity band effects, where site positions and 
the "bare" electron energies are randomly scattered, we took into account that the polaron 
shifts may be widely distributed as well. It results in a significant modification of the law 
that governs the low-temperature electrical conductivity in such materials. In the absence 
of Coulomb correlation effects, the standard Mott dependence should be substituted by 
another, non-activation one, given by Eq. @, which describes the variable-range variable- 
barrier hopping of lattice or spin polarons. 
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APPENDIX A: POLARON HOPPING RATE 

Here we will derive Eq. ( [TT]) for the polaron hopping rate. We will use a semi-classical ap- 
proach |L2| to describe lattice polaron hopping in the framework of a generalized molecular- 
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crystal model JT^]. We will further generalize this model order to take into account that the 
initial and final states may have different polaron shifts. 

In the Holstein's occurrence probability approach, the polaron hopping rate between the 
initial state i and the final state j can be written as |TB| , 



Hi 



CO 

J dEP(k 



E 



5 (tj - £i)5 (kj 



E 



(Al) 



Here 



* = t t (0) - A< n )x n, e j = E j (0) - 



(A2) 



71=1 



71=1 



are the local electron energies of the initial an final states coupled to local vibrational modes 
x n and y n with An being the coupling constants (Li is the total number of the vibration 
modes coupled to the localized carrier at Z-th site) and 



£=fe- e ^= ej (A3) 
being the time rate of change of the relative electronic energies evaluated at the event when 
the electronic terms coincide. In Eq. (|Al[ ), the brackets indicate a thermal average over 
all the configurational coordinates x n and y n and the velocities x n and V n \ P (in is the 
probability that the carrier will perform a hop at the coincidence event. 

By using the integral representation of the 5-functions and the harmonic approximation 
for the vibrational energy (see Ref. [|13| for details) one can perform averaging in Eq. (|Al|) 
and rewrite it as 



7ij 



exp {- [(A + 2Wi) 2 /4 (Wi + W 3 ) - J l3 ] /T) 



x dE PIE 



E 



exp 



- E 2 /4T (Witil + Wflfj 



(A4) 



where A = kj (0) — (0) is the electron energy difference between the carrier-free states; 
J%j = Jo exp (—Tij/a) is the hopping integral; 



(A5) 



71=1 



71=1 
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are the polaron shift and the squared average vibrational frequency of the l-th site, corre- 
spondingly. Here cui and kn are respectively the frequency and the stiffness constant of 



n-th vibrational mode that belongs to site I. 

In the so called non-adiabatic limit, when the hopping integral J 



0, the probability 



and therefore 



7, 



(n—ad) 
ij 



P(e) = 2nJfJhE« 1 



\Ji' 



71 



h V (Wi + Wj)T 



exp 



(A + 2W i: 



4(Wi + Wj)T 



(A6) 



(A7) 



which coincides with Eq.([TT]). Otherwise, in the adiabatic limit (see |T3| for details), 



Wi + W, 



exp {- [(A + 2W t ) 2 /4 (Wi + Wj) - J l3 ] /T) (A8) 



When the vibrational and coupling constants of the initial and final sites are the same, 
Wi = Wk = W, Qi = Qj = Q, the above expressions for the hopping rate coincide with these 
obtained in Ref. An expression for the hopping rate, which is similar to Eq. ([A7| ) but 

with a different pre-exponential factor, can be obtained for the bound magnetic polarons by 
using a procedure described in Ref. [|1(J. (This factor properly takes into account the vector 
nature of the order parameter (magnetization) and can be responsible for the giant positive 
hopping magnetoresistance that may be observed at low magnetic fields However, at 
high magnetic fields it coincides with that from Eq. ( |A7| ).) 
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